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Abstract 



In this article we study the solution of the Kuramoto-Sivashinsky equation on a 
bounded interval subject to a random forcing term. We show that a unique solution to 
the equation exists for all time and depends continuously on the initial data. 

>. 
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O ! 1 Introduction 



In this paper we investigate the existence and uniqueness of the solution to the Kuramoto- 
Sivashinsky (K-S) equation subject to a random forcing term. Specifically, the solution 
of 

du + {uxxxx + Uxx + uux)dt - dw = 0, (1.1) 

where w is a Q-Wiener process in a probability space (0,.F, P). The Wiener process w 
takes value in a Hilbert space to be specified later. The distributional derivative of w{t) 



d I represents an external random force. 

The usual K-S equation ((|I.1D without the dw term) has been studied as a prototypical 
example for an infinite dimensional dynamical system. It possesses a finite dimensional 



maximal attractor (|15|, 0, O, n^) and inertial manifold (pQ, 13, 1£, ^ 



Equation (|l.l| ) arises in the modeling of surface erosion via ion sputtering in amorphous 
materials The random forcing term in the model accounts for the fluctuations in the 
flux of the bombarding particles. 
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Herein we confine our attention to the case of u restricted to tlie interval / := (— /, /), subject 
to the given initial condition uq, and homogeneous Dirichlet boundary conditions, i.e. 



u{0,x) = uo{x) , -I < X < I , and u{t, -I) = u{t,l) = for t > 0. (1.2) 



We show that for any T > there exists a unique solution to ( |l.lD ,( |l.2D for < t < T, 
and establish a priori estimates for the solution. The approach we follow is similar to that 
for establishing existence and uniqueness for parabolic differential equations. Firstly we 
establish local existence (with respect to time) and then show that the solution remains 
bounded for any T > 0. For (1.1) these steps are preceeded by the introduction of a change 



of variable which enables us to consider, instead of the stochastic differential equation, a 
related deterministic equation. Local existence and uniqueness is then established via an 
application of a fixed point argument over a suitably defined space. 

The application of the fixed point theorem necessitaties expressing the nonlinear solution 
operator of the derived deterministic equation as a mapping from a space £ into itself. To 
achieve this we must show that the extensions of two operators, which arise in the analysis, 
are well defined. This effort is the major part of section ^. 

In section |^ we show that the local solution remains bounded for any T > which implies 
global existence of the solution. 

We begin our discussion by presenting in the next section several definitions and basic 
results which we use later in our analysis. 



2 Preliminaries 



As usual, we denote by LP{I), p = 1,2,... the closure of C°^{I) (the space of infinitely 
differentiable functions on /) with respect to the L'p{I) norm: 




Also, H^{I), k = l,2,..., denotes the closure of C^{I) (the space of infinitely differentiable 
functions on / which vanish at the the endponts) with respect to the Hq{I) norm: 

1 /2 

WfWHl^il) = (/ Ifl'dx + J^\ffdx + ...+ |/(^)pdx) . 
For convenience we use 

H := L^{I) and V := H^{I) . 



2 



To account for the temperal dependence we use the Banach spaces I/'iO, T; L'^{I)), with the 
associated norm: 

T / r \p/i V^"^ 



dt 



Note: The spaces L'P{0,T; Hq{I)) are defined analogously. 
A central role in the analysis below is played by the space E defined by 

E := L\0,T;L\I)) . 



We remark that this choice for E arises from the proof of lemma |3.l| and is dictated by the 
nonlinear term uUx- 

We begin by establishing the following embedding result which we combine with lemma 
]3 to establish the setting for the application of the Banach contraction mapping theorem 



(lemma 2.2 



Lemma 2.1 For any T > we have 

L°°{0,T;H)r\L'^{0,T;V) C E , (2.1) 
and there exists a constant K , independent ofT>0, such that 

\\u\\e < K (\\u\\l^^o^T;H) + II'"IIl2{o,T;V)) ,ueE. (2.2) 

Proof: We have by the Sobolev embedding theorem, (see [Q], pg. 217), that 

H^/^{D) C H^/^{D) C L\D) 
and there exists a constant Ci > such that 

lbllL4(£,) < Ci||t;||^i/2(£,) . (2.3) 

Using the interpolation inequality for H^I'^{D) in terms of L'^{D) and H^{D) we have for 
some constant C2 > and all t € [0, T] 

\\ummf^(D) < C2 hm'/'^D) Mt)\\]i^D) ■ (2-4) 

Raising both sides of ( |2.4D to the fourth power and integrating ( |2.4D over the interval [0, T], 
equation ( |2.2| ) follows using standard inequalitites and the definitions of the norms, with 
K = Ci C2/2 . 

■ 

Essential to establishing the local existence is the following contraction mapping theorem. 
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Lemma 2.2 Pg. 290) Let T denote a transformation from a Banach space £ into 

£ , d an element of £ and a > a positive number. If T (0) = 0, ||a|| < and 

\\T {zi) - T [z2)\\£ <]^\\zi - Z2\\e for WziWs < a, \\z2\\e < a, (2.5) 
then the equation 

z = d + F {z), z G f , (2.6) 
has a unique solution z ^ £ satisfying \\z\\£ < a. 



Below we use the following lemma and corollary, which describe the regularity of the solution 
to a negative self-adjoint operator. 



Lemma 2.3 f^uj] pg. 4^4) Assume that A is a negative self-adjoint operator on H and 

V = D{{-Af/^) CHCV . 
Then A and S{t) = e^^ has a continuous extension from V to V' . If 

y{t) = y{t-g) = S{t)yo + f'e^'-'^'' g{s) ds ,t G [0,T], 

Jo 

for yo G H, and g G L'^{0,T; V), then 

yeL^{0,T;H)nL\0,T;V), 
and for some constant L > 0, independent ofT>0, 

\\y\\L'^{0,T;H) + \\y\\L2{0,T;V) < -^(l|yo|k + I|5|1l2{0,T;V") ) ' (^•''') 



Corollary 2.1 For A, S{t), and yo as described in lemma 2.5, we have that 

\\S{t)yo\\E <8KT^/^ [\\S{t)yo\\U^o,T;H) + II^Wyo|lioo(o,T;V))'^' • (2- 
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Proof: Using (2.3), (|2.4|), (^) and, for notation convenience, u := S{t)yo, we have that 



r M% dt < {ciC2f r Ml 

Jo Jo 
< (C7iC2)V2 



U\\y dt 



u\\% dt + \\u\\y dt 



< (Ci C2) T/2\^ ll^llL°°(0,T;/i') + II^IIl°°(0,T; 



V) 



Taking the fourth root of both sides yields (2^) for K = Ci C2/2 



Note: From lemma 2.3 we have that S{t) yo G L°°{0, T; H) n L2(0, T; V) which guarantees 
that the right hand side of ( ^I^ ) is finite. 



3 Local Existence and Uniqueness 

Our first step in establishing local existence and uniqueness of the stochastic differential 
equation is to introduce a change of variable to reduce (1.1) to a deterministic equation. 



Denote by A the self-adjoint operator 



Au :-- 



CU . 



(3.1) 



We assume that c is chosen sufficiently large such that A is a strictly negative operator on 
the space Hq{I). 

Observe that as A is a strictly negative, self-adjoint, operator we can define (—A)" via 
Fourier analysis, with domain D({—A)°') = Hq"{I). (See |18| pg.55 for details.) 



In view of (|3.l| ) note that (|l.l| ) can be rewritten in the form 

du = {Au — uu^ + cu) dt + dw , 



(3.2) 



where the Wiener process w takes value in the separable Hilbert space H = L'^{I) and it 
has the covariance operator Q. With S{t) := e*^, t > 0, we define WA{t) via the stochastic 
integral 

WAit) := / S{t-s)dw{s). (3.3) 
Jo 

Using the substitution 



y{t,x) := u{t,x) — WAit,x) , t £ [0,T] P-a.s. , 



(3.4) 
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( |1.1| ) reduces to the deterministic problem 

yt = Ay - {y +WA){y +wa)x + c{y +wa) , (3.5) 

subject to 

2/(0, x) = uo{x) and y{t, -I) = y{t,l) = 0. (3.6) 



Note: The assumption that dw in ( |1.1| ) denotes a Q-Wiener process, together with the fact 
that ^ is a strictly negative self-adjoint operator, ensures that WAit) given by ( |3.3D has a 
version which is Holder continuous with values in D({—A)°^) for < a < 1/4, with Holder 
exponent less than (1/4 — a), (see |Q, pg. 60). Thus, with a = 1/8, in view of (|2.3|), we 
conclude that WA{t) has a continuous version in L^{I). Below we take WA{t) to denote this 
continuous version. 



The solution y satisfying ( p. 51 ) may be expressed in integral form as 
y{t) = S{t)uo + 

/ S{t- s)[-iy{s) + WAis))iyis) + ii;A(s))x + c(y(s) + wa{s))] ds (3.7) 
Jo 

= S{t)uo + F(y + WA){t) , [0,r] . (3.8) 

In the following we show the existence and uniqueness of the solution y to this integral equa- 
tion ( |3.8| ). This gives a so-called (mild) solution u for the stochastic Kuramoto-Sivashinsky 
equation (|1.1D. This is the definition of 'solution' used in this paper. 



In ( |3.8D F : E ^ E is a continuous extension of the operator 

Fo : C\[0,T]-V) ^ E 



defined by 

where 

is given by 



(Fon)(i) = f* S{t-s){Gou){s) ds ,t€[0,T], (3.9) 
Jo 

Go : C\[0,T]-V) ^ E 



(Gon)(t) = -uu^{t) + cu{t) ,t e[0,T]. (3.10) 
In view of (|3.8|) , and assuming that F is well defined — a non-trivial point whose discussion 



occupies the later part of this section — , on applying lemma p.2| we have local existence 
of the solution to (|0),(0). 



Note: The value of r for which we establish local existence and uniqueness of the solution 
on the interval [0, r], depends upon the particular realization. 
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Theorem 3.1 For uq in H there exists a random variable r taking values W-a.s. in (0,T] 
such that equations jj. j| j( |jT^j have a unique solution u on the interval [0, r]. 

Note that by a general result in |||], page 72, the solution u has a measurable modification. 
In the following the solution u refers to this measurable version. 

Proof: Observe that with z(f) = y{t) + WA{t) — S{t)uo, equation (|3.8| ) may be rewritten 
as 

z = d + T (z) (3.11) 

for d = WAit), and (z) = F{z + S{t)uo). Thus the existence and uniqueness of the 
solution to ( p. 81 ) is equivalent to that for ( |3.1ip . 

Let a = 1/6M, and ri be given by 

(6M[c(20^/' + l6Ki\\S{t)uo\\U^,^T;H) + ll^(i)^o|lioo(o,T;y))'/' 



n 



for K defined in lemma 2.1, and AI defined in lemma p 



The Ti is well-defined and it is so chosen that it will 

guarantee that .7-" is a contraction mapping (see below). 

As WA{t) is continuous with wa{0) = 0, there exists T2 such that 



^ \\wA{s)\\i4(^j)ds < a/2, for0<t<r2. 



Let r := minjri , T2} and analogous to the definition for E introduce £ as 

£ := L\0,t;L'{I)) . 



(3.12) 



With zi and Z2 satisfying < a ( = 1 / 6 M ) for i = 1,2, we have using lemma 3.2 

(12. 81), and the definition of r 



\T (zi) - {Z2)\\g < m(||zi + S{t)uo\\£ + \\Z2 + S{t)uo\\£ + c{21t)^/^ J \\Zl 

< M(||zi||^ + \\Z2\\£ + 2\\S{t)uo\\£ + c(2/t)1/4) \\zi - Z2\\g 
^2|lf • 



Finally, applying lemma |2.2| we establish the existence and uniqueness of z(t), and conse- 
quently y{t), on the interval [0,r]. 
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What remains is to establish the regularity result used for F in the proof of theore m |3.1| . 



However we must first show that F is well defined by showing Go and Fq defined by ( 3.10 ) 
and (|3.9|) have appropriate extensions. 



Lemma 3.1 The operator Gq defined by j S.lOJ can he continuously extended to 



G: E ^ L2(0,r;y') , 

satisfying 

\\G{U) - G{v)\\L2(o^T;V')<'i^^'\\\u\\E + \\v\\e + c(2/r) ^4) ||^ _ ^;||^ U^VGE. 

Proof: Let u, v, if) ^ L^(0,T;F). Denoting the duality mapping between L^(0,T;y) and 
L2(0,r;y') by (•,•), we have 

{Go{u) — Go{v),i(j) = / / {—uux +vvx)4' + c{u — v)il) dx dt 



JI 

T 



I.e. 



-{u + v){u — v)iIjx + c{u — v)'4} dx dt ^ 
Ji 2 



\{Go{u)-Go{v),ij)\ < { I j{\u\ + cf{u - vfdxdt] 



Ji 



-IV',! + \^\f dx dt] 

JI 



< I / / (\u\ + \v\ + Icl)^ dx dt] 



JI 



T 



\ 1/4 

{U - V)^ dx dt ||^||l2{0,T;V) 



Ji ) 

< (27(11^111 + \\v\\% + liclil))'^'. 

\\{U - v)\\e ||V'IIl2(0,T;V) 

< TJ^I\\u\\e + \\v\\e + c(2lTfl^) ■ 

||(n - v)\\e ||V'llL2(0,T;y) 

from which the result follows. 

For the extension of Fq and the regularity of F we have: 
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Lemma 3.2 The transformation Fq described by ( [3.^ can be continuously extended to 
F : E ^ E. Moreover there exists a constant M > 0, independent ofT>0, such that 

\\F{u) - Fiv)\\E < Mi\\u\\E + \\v\\e + c(2/r)^/4)||(n - v)\\e , u, v e E . (3.13) 



Proof : In view of the definition of Fq in ( |3.9[ ) and lemma |2.3| we have that 

F{u){t) = y{t;G{u)) e L^{0,T;H)nL\0,T;V) . 
Moreover from lemma |2.1| we have 

\\F{U) - F{v)\\e < K [\\y{-,G{u)) - y{.;G{v))Uoo^o^T;H) + 

\\y{-,G{u)) - y{-,Giv))\\ 
< i^'L|[G(n) — G(t>)||i2(o,T;F') ) using lemma |2.3| , 



< M{\\u\\e + \\v\\e + c(2/r)^/^)||(n - v)\\e , using lemma ^ , 

for M = 27^/*KL . 



4 Global Existence 



We now extend the local existence of theorem 3J. established in the previous section to 
global existence. Local existence establishes that the solution, n, lies in the solution space, 
E^ for some initial time period. We establish global existence by showing that for any time 
T the S-norm of u is finite and hence u still lies in the space E. To do this we first establish 
that the solutions are continuous with respect to the initial data. This enables us to restrict 
our attention to showing that strong solutions remain bounded in E. 



Following directly the proof of lemma 3.1 with the inner product (only) taken over the 
spatial domain, /, we have: 



Corollary 4.1 The operator G defined in lemma satisfies for u, v ^ E 

\\G{u) - G{v)\\v'<27^/''{\\u\\mr) +Mlhi) +421)'^") IW-v^Hn ' ^"^'^^ 



Lemma 4.1 The solution, y{t), of ( ^-90 depends continuous on the initial data uq G H, 
and the random forcing term WAit) G E. 
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Note that the continuous dependence on WA{t) is needed in the proof of the next lemma 
where we approximate WA{t) by regular processes. 

Proof: Let yo, and yi denote solutions of ( |3.8|) generated by uq, w'^lt), and ui, w\{t), 
respectively. Then, on the common existence interval of yo and yi, 



2/0-2/1 = S{t){uo - ui) + 



S{t-s)G{yo + w\){s) - S{t-s)G{yi+w\){s) 



ds 



From lemma ^ we have yo, and yi € L'^{0,T;V) thus (yo — yi)(0 G ^lA* o-c, (i.e. for 
almost all t E (0,T)). Using the continuity of S{t) and lemma we have that there exits 
constants Li and Ci such that 

II2/0 - 2/i||v < Li\\uo - uiWh + [ Ci\\G{yo + WA){s) - G{yi + w\)\\v' ds 

Jo 

< Li\\uo - ui\\h + 

Ci r27i/^(||yo + «;^||L4(,) +\\y,+uj\h,^j) +c{2l)'/^)\\{yo + w'i) - {yi + w\)\\^,^^^ ds , 
/i a.e. 

Using the Sobolev embedding theorem, the existence of yo, and yi € E, implies there exists 
constants C2 and C3 such that 

II2/0 - 2/iI1l4(/) < ^^2 (il^io - uiWh + \\w\ - w\\\e) + C3 \\yo - yi\\Li(i) ds 11 a.e. 
Applying Gronwall's inequality then yields 

II2/0 - 2/i||l4(/) < (ll-uo - ui\\h + \\wa - w\\\e^ e*^3t ^ Q_g_ 

from which the stated conclusion follows. ■ 
Next we establish appropriate norm estimates for the solution. 



Lemma 4.2 Let uq G H, wa be given by ( [g.^j j, and y denote the solution of 

y{t) = S{t)uQ + F{y + waW) , t G [0,r] . (4.2) 

Then, y satisfies 



sup \\y{t)\\l<\\uo\\lJo f'^^)''^ + r elo f(^)''^ gir)dT , (4.3) 
?fO,Tl Jo 



ie[o,T] 
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^ \\y{t)fydt<\\uo\\l + sup \\y{s)\\l r f{r)dT + g{T) dr , (4.4) 
se[o,T] "'0 Jo 



1 / , 3 ^ ^ 2 , ^ 11 ., ||4 



and 



where 

fit) = l(^2c + (2 + ^CiC2r + CiC2||t£;A|li4 
and 

= cWwAfn + ^II^A|li4 . 

Proof: As D{A) and C(0, T; Hq{I)) are dense in if and i?, respectively, and from lemma KTI 
we have established continuous dependence of the solution, it suffices to establish (|4.3| ),(p7^ 
for the (strong) solution of the differential equation 

^ = Ay{t) - {y{t) + WA{t)){y{t) + WA{t))^ + c{y{t) + WA{t)) , (4.5) 
2/0 = -"0 • 

We first show that 

l^umi + \mfv<f{t)\\yfH + c\\wa\\i + . (4.6) 



Multiplying (^]^) by y{t) and integrating over / we obtain 
Id 

2^\\y\\H + WVxxWh - WVxWh = - J ^y{y + WA){y + WA)xdx + cj ^ ywAdx . (4.7) 



We obtain a lower bound for llyxxlliT via: 



\\yx\\H = f yxyxdx = - [ vxxydx 
j-i J-i 



1/2 /I \ 1/2 



< ( / ^ yl^ dxj l^j ^ dxj 



- , 3r< r< \ WyxxWff + \\y\\H 



I „2 , (2 + IC1C2) 

(2 + IC1C2 

Rearranging yields 



2 



WyxxWl > (2 + Iw,) WyAl - ^ y^^^^' llylll. • (4.8) 
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Next consider the first integral on the r.h.s. of ( |4.7] ) in three pieces: 

rl rl I d 



yyyxdx = - - — {y^)dx = 0, (4.9) 
-/ J -I 6 ax 

I 1 11 

ywAWAxdx\ = |- / yxw\dx\ < -\\y\\v + -||u;A||i4(4-10) 
2, J -I I » 



and 



j ^ {yyxWA + wa^ dx = j ^yy^WAdx. 



This last term is estimated as follows. 



^yyxWAdx\ < ij^yldxj ij^y'^WAdxj 



< 



J ^y^dx] ij^y^dxj ij^w\dxj 



< lly||y||y|lL4||^A||L4 

< Ci C2\\y\\T \\y\\T Ml^ ( using {2.3),I0 ) 

< ^-^\\yfv + ^MllMlh (4.11) 

( using Young's Inequality ). 



For the remaining term on the r.h.s. of ( [1.7| ) we have 



c 



y WA dx\ < c ij y^ dxj ij w\ dx ] 



< IhWl + ■ (4.12) 



Combining (O) — (4.12) yields 



The estimate for \\y\\H in (|4.3|) now follows from the observation that (4.6) is a first order 



differential inequality for ||y||H. The bound involving ||y||y in (4.4) is established by inte- 



grating (4.6) from to T. 



We are now in a position to establish the global existence of the solution. 

Theorem 4.1 For uq £ H = L?'{I), there exists P a.s. a unique solution u{-,x) £ E of 
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Proof: From theorem p| w e have existence of the solution u{-,x) G -E, P a.s., for the 
interval [0,r]. In view of ( ^.41 ) and lemma 4^ we conclude that u{t,x) remains bounded in 
E, P a.s. for all t > 0, which implies global existence of the solution to (1.1), (|L^). ■ 



Finally we remark that by following the same argument as in Brannan et al. we can 
show that the solution is actually Holder continuous in space with exponent less than g. It 
is also possible to consider multiplicative noise in equation (1.1). The approach in this paper 
should also apply to other similar parabolic type stochastic partial differential equations. 



References 

[1] R. Adams, Sobolev Spaces, Academic Press, New York, 1975. 

[2] J. Brannan, J. Duan and T. Wanner, Dissipative Quasigeostrophic Dynamics under 
Random Forcing, J. Math. Anal. Appl. 228 (1998), 221-233. 

[3] P. Collet, J.-P.Eckmann, H. Epstein and J. Stubbe, "vl global attracting set for the 
Kuramoto-Sivashinsky equatiori\ Comm. Math. Phys. 152 (1993), 203-214. 

[4] R. Cuerno, H. A. Makse, S. Tomassone, S. T. Harrington and H. E. Stanley, ^^Stochastic- 
model for surface erosion via ion sputtering: dynamical evolution from ripple morphol- 
ogy to rough morphology,''^ Physical Review Letters, 75, (1995) 4464-4467. 

[5] G. Da Prato and J. Zabczyk, Stochastic Equations in Infinite Dimensions, Cambridge 
University Press, 1992. 

[6] G. Da Prato and J. Zabczyk, Ergodicity for Infinite Dimensional Systems, Cambridge 
University Press, 1996. 

[7] J. Duan, H. V. Ly and E. S. Titi, "T/ie effect of dispersion on the dynamics of the 
Kuramoto-Sivashinsky equation," preprint, 1995. 

[8] J. Duan and V. J. Ervin, "Dynamics of a Nonlocal Kuramoto-Sivashinsky Equation," 
J. Diff. Eqns., 2, (1998), 243-266. 

[9] N. M. Ercolani, D. W. McLaughlin and H. Poitner, "Attractors and transients for a 
perturbed periodic KdV equation: a nonlinear spectral analysis," J. Nonlinear Sci., 3, 
(1993), 477-539. 

[10] C. Foias, B. Nicolaenko, G. R. Sell and R. Temam, "Inertial manifolds for the 
Kuramoto-Sivashinsky equation and an estimate of their lowest dimension," J. Math. 
Pures et Appl. 67 (1988), 197-226. 



13 



[11] J. Goodman, ^^Stability of the Kuramoto-Sivashinsky and related systems," Comm. Pure 
Appl. Math. 47 (1994), 293-306. 

[12] J. S. Il'yashenko, "Global analysis of the phase portrait for the Kuramoto-Sivashinsky 
equation," J. Dynamics Diff. Eqn. 4 (1992), 585-615. 

[13] M. S. Jolly, I. G. Kevrekidis and E. S. Titi, "Approximate inertial manifolds for the 
Kuramoto-Sivashinsky equation: analysis and computations," Physica D 44 (1990), 
38-60. 

[14] K. B. Lauritsen, R. Cuerno and H. A. Makse, "Noisy Kuramoto-Sivashinsky equation 
for an erosion model," Physical Review E, 54, (1996), 3577-3580. 

[15] B. Nicolaenko, B. Scheurer and R. Temam, "Some global dynamical properties of a 
class of pattern formation equations", Comm. in PDEs, 14, (1989), 245-297. 

[16] A. Pazy, Semigroups of Linear Operators and Applications to Partial Differential Equa- 
tions, Springer- Verlag, 1983. 

[17] J. C. Robinson, "Inertial manifolds for the Kuramoto-Sivashinsky equation," Phys. 
Lett. A 184 (1994), 190-193. 

[18] R. Temam, Infinite- Dimensional Dynamical Systems in Mechanics and Physics, 
Springer- Verlag, New York, 1988. 

[19] R. Temam and X. Wang, "Estimates on the lowest dimension of inertial manifolds for 
the Kuramoto-Sivashinsky equation in the general case," Diff. Integral Eqns 7 (1994), 
1095-1108. 

[20] E. Zeidlcr, Nonlinear Functional Analysis and its Applications II/ A, Springer- Verlag, 
New York, 1990. 



14 



